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Abstract. Near IR spectra at low (R~800) and medium (R~9000) resolution, obtained with ISAAC at VLT, have 
been used to pose constraints on the evolutionary state and accretion properties of a sample of five embedded 
YSOs located in the R CrA core. This sample includes three Class I sources (HHIOO IR, IRS2 and IRS5), and two 
sources with NIR excesses (IRS6 and IRS3). IRS5 and IRS6 have been discovered to be binaries with a separation 
between the two components of 97 and 78 AU, respectively. Absorption lines, typical of late-type photospheres, 
have been detected in the medium resolution spectra of all the observed targets, including HHIOO IR and IRS2 
which have high values of infrared continuum veiling (vk = 6 and 3, respectively). These two sources also present 
low resolution spectra rich of emission lines (HI, CO and plenty of other permitted lines from neutral atoms) likely 
originating in the disk-star-wind connected regions. Among the features observed in HHIOO IR and IRS2, Na I 
at 2.205 fim and CO at 2.3/^m, which are more commonly used for stellar classification, are detected in emission 
instead of absorption. Several strong photospheric lines, which lie around 2.12 and 2.23fim and whose ratio is 
sensitive to both effective temperature and gravity, are proposed as independent diagnostic tools for this type of 
sources. We derived spectral types, veiling and stellar luminosity of the five observed sources, which in turn have 
been used to infer their mass (ranging between 0.3-1.2 M0) and age (between 10^ and 10® yr) adopting pre-main 
sequence evolutionary tracks. We find that in HHIOO IR and IRS2 most of the bolometric luminosity is due to 
accretion (Lacc/Lboi ~0.8 and 0.6 respectively), while the other three investigated sources , including the Class 
I object IRS5a, present a low accretion activity {Lace/ Lt,oi < 0.2). Mass accretion rates of the order of 2 10~® 
and 3 10~^ M© yr~^ are derived for HHIOO IR and IRS2, respectively, i.e. higher by an order of magnitude with 
respect to the average values derived in T Tauri stars. We observe a general correlation between the accretion 
luminosity, the IR veiling and the emission line activity of the sources. In particular, we find that the correlation 
between Lace and Lsr-y, previously reported for optical T Tauri stars, can be extended also to the embedded 
sources, up to at least one order of magnitude larger line luminosity. A correlation between the accretion activity 
and the spectral energy distribution slope is recognizable but with the notable exception of IRS5a. Our analysis 
therefore shows how the definition of the evolutionary stage of deeply embedded YSOs by means of IR colors 
needs to be more carefully refined. 
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Line: formation 

velopes. In T Tauri systems, where the contribution from 
the dusty envelope has become negligible, there is the pos- 
sibility to disentangle, through observations at different 
wavelengths, the emission from the star itself from that 
of the different active circumstellar regions, and to define 
in detail the properties of these latters. Hence the stellar 
parameters of T Tauri stars have been derived through 
the observations of their optical and IR photospheric fea- 
tures (e.g.Luhman & Rieke 1998), their accretion region 
revealed through the observation of the UV excesses over 
the stellar photosphere (e.g. GuUbring et al. 1998), and 
finally their coUimated jets probed down to close to their 
base through high resolution observations of optical for- 



1. Introduction 

During the first evolutionary stages of a young star, the 
protostellar system is characterized by a variety of phe- 
nomena occurring in different regions of the circumstel- 
lar environment. These include the accretion disk, the 
magnetospheric accretion region, coUimated jets and, in 
more embedded sources, the still not dispersed dusty en- 
Send offprint requests to: Brunella Nisini, 
msini@mporzio.astro.it 

* Based on observations collected at the European Southern 
Observatory, Chile (ESO N.069.C-0269, 070.C.-0138, 063.F 
0691) 
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bidden lines (e.g. Bacciotti et al. 2000). The possibility to 

separate the different emission components contributing 
to the stellar system has greatly advanced the comparison 
of observations with models describing the accretion and 
ejection mechanisms in these stars, as well as the disk evo- 
lution and dispersal. It is however fundamental to under- 
stand how and if these models works for less evolved sys- 
tems, where the circumstellar activity is much higher and 
the accretion luminosity is expected to dominate over the 
luminosity of the still forming star. To this aim, a sample 
of interesting objects to bo investigated are sources with 
no optical counterpart, characterized by steeply rising IR 
SEDs and with infrared colors indicating the presence of 
intrinsic excesses due to circumstellar dust (the so-called 
Class I sources). The standard paradigm is that most of 
these sources are protostars still embedded in their origi- 
nal infalling envelope and expected to derive most of their 
luminosity from accretion through a massive circumstellar 
disk. For such more embedded sources, however, the large 
extinction coupled with strong emission excesses due to 
the circumstellar envelope, make the investigation of phe- 
nomena occurring very close to the protostar extremely 
difficult and thus their real nature remains often obscure. 
Infrared spectroscopy in this respect can significantly im- 
prove the study of the close environments of these sources, 
since different gas features which arc present in the near 
IR atmospheric windows can be related to different emit- 
ting regions. The near IR spectra of these objects are often 
characterized by the presence of ionized (e.g. Fc+) , neu- 
tral (e.g. HI, He) and molecular (CO, H2) emission lines. 
Such lines originate in the inner and warm regions of the 
disk, in the accretion region connecting the disk to the 
star, and in the related phenomenon of winds and jets, and 
consequently they have been used to probe such different 
environments (e.g. Nisini et al. 2002, Najita et al. 1996, 
Davis et al. 2002, Muzerolle et al. 1998). The large values 
of the IR excess characteristic of these source, however, 
produce a large veiling which makes extremely difficult 
to infer their stellar properties by means of their under- 
lying photospheric features and therefore to correlate the 
spectral probes of circumstellar activity with the physical 
parameters of the star itself. 

Recently, the first direct observation of absorption lines 
from the photosphere of a Class I source in the p Oph cloud 
has been reported by Greene & Lada (2002). The observed 
features are typical of late type stars, confirming the low 
mass nature of these sources, and have been nscd to deter- 
mine for the first time the stellar properties of one of such 
highly obscured objects, which resulted to have ^ 70% of 
its luminosity due to accretion. The potentiality offered by 
high-resolution near IR spectroscopy to reveal the central 
protostellar source characteristics opens in turn the pos- 
sibility to both derive the mass and ages of these young 
objects, as well as infer quantitatively the mass accretion 
rate and accretion luminosity to be compared with those 
of expected more evolved soiirces. In order to make this 
comparison significant it is also important to have this 
information in a sample as large as possible of embed- 



ded sources, in such a way to better constrain their global 
accretion properties. In this paper we report the results 
of an analysis of medium/high resolution near IR spec- 
tra of a sample of embedded young stellar objects of the 
nearby R CrA star forming cloud (D=130 pc, Marraco & 
Rydgren 1981). The selection of sources from the same 
cloud assures that the properties of the sources reflect in- 
trinsic differences and not the different characteristics of 
the hosting cloud. The aim of this paper is twofold: 1) to 
infer, through the observations of photospheric features, 
the stellar physical parameters of the considered objects 
and from these indirectly derive their age, accretion lumi- 
nosity and mass accretion rates. This information will be 
used to understand to which extent the spectral energy 
distribution of embedded IR sources is really indicative 
of their evolutionary stage; 2) to establish if there is a 
correlation between the accretion luminosity of the em- 
bedded stars and the occurrence and luminosity of their 
IR emission lines in order to have an indirect instrument 
to quantitatively derive the accretion rate from simple ob- 
servational tools. We describe the considered sample in §2 
and the observations obtained at VLT with the ISAAC 
instrument in §3. In §4 and §5 the characteristics of the 
emission line spectra are described, while in §6 and 7 we 
present the procedure adopted to derive the stellar and 
circumstellar properties from the detected absorption fea- 
tures. Finally, in §8 we discuss the implications of our 
findings. 



2. Description of the selected sample 

The soiirces investigated in this paper have been taken 
from the sample of IR young stellar objects in the R CrA 
cloud core listed by Wilking et al. (1997). Prom this list, 
we have selected all the IR sources with K magnitudes 
smaller than 10.5 mag and IR colors compatible with the 
locus of young embedded protostars with intrinsic IR ex- 
cesses (Lada & Adams 1992). The list of sources is pre- 
sented in Table 1 with their coordinates, IR magnitudes 
and luminosities, when available. All the sources have lu- 
minosities ranging between 1 and 15 L©, and thus they 
are expected to be low mass young stellar objects. 
Three of the considered sources, namely HHIOO IR, IRS2 
and IRS5, can be classified as Class I protostars on the 
basis of their spectral energy distribution (having spec- 
tral index between 2 and lO/zm > 0, Wilking et al. 1986). 
These three sources are known to have large dusty en- 
velopes as testified by the presence of strong ice bands and 
silicate features in their 3, 5 and lO/zm spectra (Brooke et 
al. 1999, Whittet et al. 1996, Alexander et al. 2003, Chen 
& Graham 1993, Pontoppidan et al. 2003). HHIOO IR has 
been also suggested as the driving source of the Herbig 
Haro objects HHlOl and HH99 (Hartigan & Graham 
1987) and of a bipolar molecular outflow (Anglada et al. 
1989). 
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Table 1. Observed sample 



Source 


a [ZUUU) 


o [Z\j\j\j) 


A 


(K-L) {H-K) 
magnitudes 


[J-M) 




xiei. 


T 

(T ^\ 


IRS2 (TS13.1) 


19 01 41.5 


-36 58 29 


7.16 


2.72 


2.73 


4.07 


0.98 


TS 


12" 








7.2 




2.7 


3.9 




SOFI 




IRS5 (TS2.4) 


19 01 48.0 


-36 57 19 


10.33 


2.75 


3.32 


3.89 


1.12 


TS 


3" 


IRS 5a 






10.9 




3.8 






ISAAC 


2 


IRS5b 






11.5 




3.4 






ISAAC 




IRS6 (TS2.3) 


19 01 50.6 


-36 56 38 


10.34 


1.74 


2.37 


3.64 


~0 


TS 


-0.6'' 


IRS6a 






11.1 




2.2 






ISAAC 


~0.3 


IRS6b 






11.7 




3.1 






ISAAC 




HHIOO IR (TS2.6) 


19 01 50.7 


-36 58 10 


8.04 


3.36 


3.5 


4.87 


1.47 


TS 


14" 








7.4 




2.5 


4.1 




SOFI 




IRS3 (TS4.4) 


19 02 05.1 


-36 58 56 


8.64 




2.02 


2.26 




TS 


-0.3"^ 








9.50 




0.93 


0.14 




B92 





"Spectral index {dlog{X F\) /dlogX) between 2 and 10/xm. 

''Origin of photometric data: TS = Taylor & Storey (1984), SOFI = our NTT-SOFI observations (March 2003), ISAAC = 
VLT-ISAAC archive data (1999), B92 = Burton et al. (1992). 
'^From Wilking et al. 1992 

''estimated from the photometry up to the Q band given by Wilking et al. 1986 



3. Observations 

Observations were carried out at VLT/UTl with the 
ISAAC spectrograph in SW mode on 12 and 13 July, 2002. 
At the spatial resolution provided by ISAAC (i.e. 0.14 arc- 
sec/pixel), the sources IRS5 and IRS6 were recognized to 
be binaries, with a separation of ~0.75 and 0.6 arcsec (i.e. 
97 and 78 AU), respectively. Figure ^ shows H and K 
images obtained from VLT pre-imaging of these stars, in 
which the two components that we named a and b are 
clearly visible. Due to their alignment in the sky, a single 
slit, with a position angle of 32.6°, was adopted to observe 
simultaneously IRS5a and b, and IRS6a. 

Low resolution spectra through the J(l.l-1.4/im), SH 
(1.4-1.82/im) and SK (1.82-2.5/xm) filters were acquired 
with a slit aperture of 0'.'6 corresponding to a resolution 
from 750 (in the K band) to 960 (in the J band). For each 
source, medium resolution spectra were also acquired in 
two spectral segments centered at 2.261yLtm and 2.161/.tm 
each covering — 0.122/im around the central wavelength. 
HHIOO IR, IRS2 and IRS3 were also observed in an 
additional segment in the H band, centered at 1.6295/xm 
and covering '^0.080/im. For these observations we used 
the 0'.'3 slit correspondig to a nominal resolution of 
R~8000 and R~9000 for the K and H observations, 
respectively. Spectra of standard stars of O spectral type 
were also obtained in both low and medium resolution 
at an airmass similar to the scientific spectra to cor- 
rect for telluric absorption. The telluric spectra were 
carefully cleaned by any intrinsic HI and He absorption 
feature before being used. In the low resolution spectra, 
wavelength calibrations were performed using a Xenon 
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24.0 
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Fig. 1. VLT-ISAAC H and K images of the sources IRS5 
and IRS6 were the discovered binary components a and b 
are indicated. 



lamp spectrum taken at the end of the night. In medium 
resolution, the lamp calibration was refined each time on 
the OH sky lines observed in the spectra. In addition, the 
object and standard spectra were carefully aligned each 
other against the atmospheric absorption features before 
being rationed. This procedure leads to a wavelength 
calibration error in medium resolution mode that we 
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estimated of the order of 0.1 A(i.e. about 1-2 kms ^). 

The spectra were obtained during not photometric 
nights, with a seeing variable between 1.5 and 2 arcsec. 
To flux calibrate our spectra we therefore used the broad 
band photometry obtained at SOFI with the NTT instru- 
ment during March 2003, whose values are reported in 
Table 1. For the IRS5 and IRS6 sources, since the NTT 
resolution is not high enough to resolve the binaries, we 
derived the magnitudes of the relative components using 
H and K VLT-ISAAC images taken from the ESO archive 
(ID programme 63.1-0691, April 1999). These are also re- 
ported in Table 1. From this table, and from the compari- 
son of our derived photometry with the photometry given 
by Taylor & Storey (1984), we notice that a certain degree 
of variability is present in the considered sources, with the 
largest variation exhibited by HHIOO IR and IRS3. HHIOO 
IR source was already known to be highly variable, with 
near IR magnitudes changing by values as large as more 
than two magnitudes on timescales of one year (Molinari 
et al. 1993). 

The complete low resolution spectra of the targetted 
sources are shown in Fig. 13 The NIR SEDs of all the 
sources but IRS3 are characterized by a continuum steeply 
rising from the J to the K band. The IRS 3 SED peaks in 
the H band and then declines. This evidence, and the fact 
that the H-K color of this source more recently measured, 
is smaller by about 1 mag with respect to the value given 
by Taylor & Storey (1984) (see Table 1), indicate that the 
IRS3 IR excess is significantly reduced with time. In the 
low resolution spectra of IRS3, IRS5a,b and IRS6a some 
features in absorption are present, while no significant line 
in emission is detected. Conversely, IRS2 and HHIOO IR 
are characterized by rich spectra in emission with no ab- 
sorption features (see also Figs. I3I4|I . Moreover in both 
of these spectra, and to a lesser extent also in IRS6a, we 
detect two broad features at ^^1.10 and 1.17/j,m which re- 
main unidentified. 

Figs. I5I6I and [3 show the normalized, high resolu- 
tion spectral segments observed in each source. Given the 
weakness of IRS5b with respect to the other component, 
it was not possible to extract its medium resolution spec- 
trum and therefore only the MR spectrum of IRS5a has 
been analyzed. All these spectra are characterized by the 
presence of many features in absorption. This is true also 
for HHIOO IR and IRS 2 where the higher resolution has al- 
lowed to detect narrow and weak features which remained 
undetected in the low resolution spectra. We will sepa- 
rately discuss the emission and absorption spectra of the 
observed sources in the following §4 and 5. 

4. The emission lines 

HHIOO IR and IRS2 arc the only sources showing a promi- 
nent emission spectrum rich of features (Figs. 3 and 4). 
Such features mainly include permitted transitions of neu- 
tral atoms. The HI lines are the strongest of them: we de- 
tect Br7, Pa/3 and other lines from the Brackett series up 
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Fig. 2. Low resolution J,H and K spectra of the observed 
sources 

to n=20 in HHIOO IR and IRS2, while a weak Br7 emission 
is also mesured in the MR spectrum of IRS5a. The HI lines 
observed in HHIOO IR have been already presented and 
discussed in Nisini et al. (2004); in Table 2 we report the 
fluxes and the kinematical properties of the HI lines de- 
tected both in the low and high resolution spectra of IRS2 
and IRS5a. A full account on the interpretation of the HI 
line emission in HHIOO IR has been given in Nisini et al. 
(2004). The Brackett decrement observed on HHIOO IR is 
consistent with emission from moderately optically thick 
lines originated from a compact region of few stellar radii. 
Moreover, the fact that the Brackett line profiles observed 
in high resolution are almost symmetrical with linewidths 
decreasing with the upper quantum number, is an indica- 
tion in favour of an origin in an expanding stellar wind 
rather than in an accretion region, as often suggested for 
these lines (e.g. Muzerolle et al. 1998). In such an hypoth- 
esis, a mass loss rate of ~10^® Mq yr^^ has been derived 
for HHIOO IR by using a simple model of an ionized wind. 
The Brackett lines observed in IRS2 shows characteristics 
similar to the HI lines of HHIOO IR, such as the symme- 
try of the Br7 line and large linewidths. The Brl2 and 
13 profiles are also narrower than the Br7 line, but they 
appear less symmetric partially due to a superposition of 
few photospheric absorption lines which make difficult to 
reconstruct the intrinsic line profile. The IRS2 Brackett 
decrement resulting after having dereddened the observed 
lines with an extinction value Ay— 20 (see Table 7) is com- 
patible with line optical depths smaller than in the HHIOO 
IR case (see e.g. Fig. 2 of Nisini et al. (2004)), indicating 
lower densities of ionized material. If interpreted in the 
same wind model adopted for HHIOO IR, such a decre- 
ment is consistent with a mass loss rate of ionized gas of 
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~10~^ M0yr~^, corresponding to a total mass loss rate 
of ^ 1-2 10^^ Mq yr~^ assuming a ionization fraction be- 
tween 0.05 and 0.1. 

The other emission lines detected in the low resolution 
spectra of HHIOO IR and IRS2 are listed in Table 3 where 
we have included only the features having an S/N larger 
than 3 and linewidth comparable to the instrumental res- 
olution. In this table we indicate, when possible, the most 
probable identification done using the NIST database^ 
and guided by the spectra of other YSOs (Kelly et al. 
1994) or objects with likely similar excitation conditions 
(Walmsely et al. 2000, Hamann et al. 1994). We tried as 
much as possible to check for the plausibility of a line iden- 
tification by verifying that all the transitions of the same 
atom expected to originate from the same term were also 
detected. The observed atomic lines include several CI, 
01, Sil and Mg I transitions, all coming from levels with 
energies < 12.3 eV. The only ionic lines detected are from 
singly ionized Fe and Ca, which have ionization potential 
lower than 8 eV. In the high resolution spectral segments, 
in addition to the Br7, Brl2 and Brl3 lines, the Na I dou- 
blet lines at 2.206 and 2.209^m, and the 2.12^m H2 v=l-0 
S(l) lines are detected in emission in HHIOO IR and IRS2 
and their fluxes and kinematical characteristics reported 
in Tab 3. 

Most of the observed atomic lines can be excited by 
both recombination and photo-excitation. For the two 01 
lines at 1.1287,1.3165 fim fluorescence excitation by ei- 
ther an UV continuum or resonant absorption of HI Ly/3 
photons have been proposed as the most likely excita- 
tion mechanisms (see e.g. Kelly et al. 1994, Walmsley 
2000). The several CI lines are on the other hand the 
brightest permitted lines predicted by recombination the- 
ory (Escalante & Victor 1990). The several lines from Sil 
and Mgl, which have low ionization potentials, may also 
originate from recombination cascade. Some of the atomic 
species observed are those responsible for strong permitted 
optical lines in T Tauri stars, such as Nal and Call. The 
levels from which the IR lines originate lie at excitation 
energies larger than the levels of the optical lines, but it is 
conciveable that the origin of the emission must be simi- 
lar. Among different possibilities, emission from the disk- 
star accretion region (Basri & Bertout 1989, Hartmann 
1998), from a stellar wind (Hartmann et al. 1990) and 
from an active cromosphere (Calvet et al. 1984) are the 
most commonly proposed. The kinematical characteristics 
of the two Na I lines detected in high resolution are re- 
markably similar to those of the Br7 emission, both in 
terms of linewidths and Vlsr, which suggests that they 
come from the same emission region. We also note the ab- 
sence of forbidden lines, e.g. from [Fell], which indicate 
that the emission line region is characterized by high den- 
sities where forbidden lines with critical densities lower 
than 10^ cm~^ are quenched. 

The weak H2 2.12^m hue detected in HHIOO IR and 
IRS2 (Table 4) may be an indicator of both shocked gas 

^ [http://physics.nist.gov/PhysRefData/ ASDl/l 



in the interaction of the stellar wind with the ambient 
medium, and emission from the inner warm region of a 
circumstellar disk. The line is blueshifted and only barely 
resolved, with an intrinsic width of ^ 20 kms^^ . Given 
these characteristics, the detected H2 emission may trace 
the base of a compact molecular outflow, such as those 
already detected in other Class I sources by Davis et al. 
(2001). In our case the emission is not spatially resolved, 
however the adopted slit may not be aligned with the 
compact jet emission. The origin in a circumstellar disk 
is however equally plausible. The upper limit on the line 
widths, together with the source mass estimates given in 
Tab 7, indicate that for a keplerian disk, the observed line 
are consistent with an emission region at > 1 AU from the 
stellar source. Ro- vibrational H2 emission with linewidths 
similar to those observed by us have been detected from 
disks of some T Tauri stars (Bary et al. 2003). In conclu- 
sion both the jet and the disk hypothesis seem to be plau- 
sible for the 2.12^m line origin. Finally CO in emission 
from different vibrational band heads in the 2.3/im region 
is also detected in HHIOO IR and IRS2. CO band head 
emission is a tracer for dense molecular gas at T~2000- 
3000 K, and it is usually interpreted as a sign for the 
presence of active accretion disks. The reason is that to 
have this band in emission large columns of warm gas need 
to be present, and accretion may represent a viable disk 
heating mechanism (Najita et al. 1996). 



The spectra of HHIOO IR, IRS2 and IRS6a show, in 
addition to the emission lines, also two broad and strong 
unidentified features in emission, at 1.098 and 1.173 /zm, 
whose parameters are given in Table 5. Such features seem 
not originated by spurious or instrumental effects since 
they are not observed neither in the IRS3 spectrum nor in 
the spectra of the standard stars. They do not appear spa- 
tially resolved along the slit, being confined on the source. 
The large widths and strengths of these features resam- 
ble the characteristic shape of PAH emission, which how- 
ever are so far never been observed at these wavelengths. 
The observations of PAH at such short wavelengths would 
imply electronic excitation of these macro-molecules and 
consequently the presence of strong UV fields which seem 
not to be associated with the low luminosity sources we 
are considering. The attribution of these features to PAH 
emission is also not favoured by the non-detection of the 
vibrationally excited PAH feature at 3.3/^m in the spec- 
trum of HHIOO IR (Whittet et al. 1996). Diffuse inter- 
stellar bands (DIBs) have been detected in the consid- 
ered spectral range (Jobhn et al. 1990) but they are fea- 
tures usually observed in absorption against a strong con- 
tinuum. Unidentified dust emission features at 1.15 and 
l.bfim have been finally observed towards the refiection 
nebula NGC7023 (Gordon et al. 2000). Such features are 
however much broader (^ 0.2-0.3/im) and more diffuse 
than those detected by us. 
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Table 2. Emission lines observed in HHIOO IR and IRS2 (excluding HI lines) 



Aobs 


F(AF) 


i^(AF) 
f s"""^ cm~^ 


Identification {Xvac in /im) 


lira 


10"^'' erg 






T TT T 1 C\C\ TFl 

HHIOO IK 


IRS 2 




1.1291 


1.52(0.05) 


0.89(0.07) 


OI 3p^P-3(fD° 1.1287 


1.1327 


0.05(0.05) 




CI 3d^D°-3p^P 1.1333 


1.1339 


0.15(0.05) 




CI 3d^F°-3p^D 1.1339 


1.1641 


0.34(0.08) 




CI 3d^D°-3p^D 1.1632, 1.1651 


1.1754 


0.59(0.08) 




CI 3d^F''-3p^Dl. 1751, 1.1756, 1.1758 


1.1831 


0.37(0.08) 




Mgl As^S-4p^P° 1.1831 


1.1839 


0.55(0.08) 




? 


1.1835 


0.84(0.08) 




Call 5p^P°-5s^S 1.1842 


1.1873 


0.27(0.08) 




? 


1.1888 


0.84(0.08) 


0.27(0.06) 


CI 4s^P°-3p^D1.1883 


1.1954 


0.32(0.05) 




Call 5/P°-5s^5' 1.1953 


1.1978 


1.02(0.05) 


0.83(0.06) 


? 


1.1992 


0.81(0.05) 


0.20(0.06) 


Sil 4p^D-4s^P° 1.9948 


1.2035 


0.60(0.05) 


0.35(0.06) 


Sil 4p^D-As^P° 1.2035 


1.2086 


0.19(0.05) 




Mgl if^F'^-M^D 1.2087 


1.2111 


0.62(0.05) 




Sil 4p^D-4:S^P° 1.2107 


1.2275 


0.30(0.04) 




Sil 4p^D-4:S^P° 1.2274 


1.2542 


0.26(0.04) 






1.2568 


0.36(0.04) 


0.3(0.1) 


CI 3d^P''-3p^P1.2566 + 
[Fell] a^Dg/2-a'^Dy2 1-2567 


1.2611 


1.93(0.04) 


1.6(0.1) 


CI 3d^P''-3p^P1.2617 


1.2640 


0.93(0.04) 


0.8(0.1) 


Fell z^Py/a-c^Ps/a 1.2646 


1.2675 


0.60(0.04) 


0.5(0.1) 


? 


1.2701 


1.02(0.04) 


1.0(0.1) 


? 


1.2784 


1.18(0.04) 


1.2(0.2) 


? 


1.3176 


0.64(0.12) 




01 3p='P-4s^5'° 1.3166 


1.4884 


10.3(1.3) 




Mgl 4fF°-2,cPD 1.4882 


1.5047 


20.1 (0.7) 


3.5 (0.9) 


Mgl 4p^P°-4s^S' 1.5044 + Br23 


1.7113 


11.1(1.8) 




Mgl Ap^P^-As^S 1.7113 


2.294 


540(30) 


87(17) 


CO v=2-0 


2.324 


600(30) 


110(25) 


CO v=3-l 


2.354 


500(30) 




CO v=4-2 


2.324 


390(40) 




CO v=5-3 


2.324 


360(40) 




CO v=6-4 





5. The absorption features 

In the two K band spectral segments, the main absorption 
features observed are those of lines from neutral atoms like 
Na, Al, Fe, Mg, Ti and Si, as well as the CO 2-0 band; 
these are the typical features expected from the photo- 
sphere of late type (K-M) stars (e.g. Wallace & Hinkle 
1996, 1997) and commonly observed in pre-main sequence 
sources (e.g. Luhman & Rieke 1998). In HHIOO IR and 
IRS2 such features are detected at only 1-5% of the contin- 
uum level. In the H band spectral segment, several more 
features are detected, most of which are very weak. The 
most significant of these features are indicated in Figs. l5l6l 
andQand have been identified through the spectra of stan- 
dard stars of Meyer et al. (1998) and the Kurucz atomic 



spectral line database ^ . The absorption lines are resolved, 
with widths of several tens kms~^ which vary depend- 
ing on the considered feature in the same objects. Such 
a broadening is due to different contributions, including 
rotation and Zeeman broadening, since these young stars 
are expected to be fast rotators and to have strong stel- 
lar magnetic fields. More specifically, the width of lines 
selected for having small Lande factors and not blended 
with other features are of the order of ~20kms~^ in IRS3, 
40kms-i inIRS2, HHIOOIR and IRS6a, and 60km s-i in 
IRS5a. This indicates that all the sources but IRS3 pos- 
sess large v sin i values, comparable with the projected 



^ 1995 Atomic Line Data (R.L. Kurucz and B. Bell) 
Kurucz CD-ROM No. 23. Cambridge, Mass.: Smithsonian 
Astrophysical Observatory 
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Table 4. Other emission lines in the high resolution spectra 
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HHIOO IR IRS2 
A Line Vlsr F±AF Vlsr A 1/ F±AF 

fim kms"'^ kms^^ 10~^^ erg s~^ cm~^ kms~^ kms~^ 10~^^ ergs~^ cm 



2.1218 Ha -6 ^20 5.4 ± 0.3 -6 <,20 3.4 ± 1.0 

2.2062 Na I ^Pa/a-^Si/a 5 142 32.0 ± 0.8 ... ... 8.8 ± 1.5 

2.2090 Na I ^-Pi/a-^gya 2 133 42.0 ± 1.5 -11 117 19.0 ± 0.5 



Table 5. Parameters of the l.l^m unidentified features 



HHIOO IR IRS2 
A F±AF FWHM F ± AF FWHM 

10"^** ergs~^ cm~^ fim 10"^* ergs~^ cm~^ 



1.098 42.9±0.9 7.410"^ 5.1±0.2 7.510"^ 

1.173 13.2±0.7 1.310"^ 1.4±0.1 1.2 10"^ 




rotational velocities estimated in other flat-spectrum and 
Class I sources and higher than the average v sin i val- 
ues measured in T Tauri stars (< vsini > ^15kms^^ , 
Bouvier et al. 1986). 

5.1. Spectral type and veiling measurement 

In the assumption that the observed absorption lines 
originate from stellar photosphere, their equivalent width 



(EW) should be sensitive to the stellar effective tempera- 
ture and gravity. In the sources we are considering, how- 
ever, the line EWs are strongly affected by the presence of 
continuum excesses, which can be quantified through the 
veiling parameter, defined as the flux at a given wave- 
length in excess with respect to the photospheric flux 
(rx = Fex,\/ F*,\). The veiling is actually a parameter 
that we also want to derive from the data, giving us in- 
formation on the circumstellar activity of the stars. The 
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1.05 7 
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X(^im) 



2.28 



2.3 



Fig. 5. Continuum normalized, medium resolution spectra of HHIOO IR. The most important emission and absorption 
features are labelled. 



ratios of different lines close in wavelengths is not largely 
affected by the veiling and can thus be used for defining 
the spectral type. Then, the absolute value of the line EWs 
can be in turn compared with the intrinsic EW of spectral 
standard lines to estimate the IR veiling. Several diagnos- 
tic lines, mainly located in the K band, have been identi- 
fied by different authors as particulary suited for spectral 
classification (e.g. Luhman & Rieke 1998). Among them, 
different features lying in the region around the Na I dou- 
blet at 2.206/xm are good diagnostic tools for the effective 
temperature, while the CO 2-0 bandhead is particularly 



sensitive to the gravity (e.g. Doppmann & Jaffe 2003). 
However, we have seen that in HHIOO IR and IRS2 both 
the Na I doublet lines and the CO band-heads are seen 
in emission instead of absorption, and therefore cannot 
be used for spectral classification. Therefore, to derive the 
stellar properties of our sources, we have used different 
and close in wavelength diagnostic lines, mainly located 
in the K band. Details on the considered lines are given 
in Appendix A. In practice to derive the spectral type and 
luminosity class, we have compared the ratios of selected 
absorption lines observed in our sources, with the same 




1 1 1 1 


1 1 1 1 1 1 1 
IRS5a 


1 1 1 1 1 1 1 

Bit 






- Jl 


Si 


TiSi Sil 


1 Na 
Na 




■ , , , 1 


. . . 1 . . . 


Ti 

. . . 1 . . . 





2.12 2.14 2.16 2.18 2.2 2.22 




\{fj.m) 



Fig. 7. As in Fig. but for IRS5a (left) and IRS6a (right) 

ratios as measured in a grid of standard star spectra, se- 
lecting the standard which minimizes the scatter among 
the different ratios. 

A difficulty in this procedure is given by the limited 
number of available spectra of standard stars at a reso- 
lution similar to that adopted for our observations. For 
our analysis we have mainly used the Wallace & Hinkle 
(1997) (for K band), and Meyer et al. (1998) (for the H 
band) standard samples, which are however at resolution 




3000, thus not always allowing to compare lines resolved 
in our spectra but blended in the standard stars spectra. 
Moreover, the Wallace & Hinkle (1997) sample lacks spec- 
tra of dwarfs stars between K5 and M2. We therefore also 
used the dwarfs standard star spectra obtained by Greene 
& Lada (2002) which have a finer grid between K4V and 
M2V at a resolution of about 15000, but do not enterely 
overlap with our observed spectral range. 
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Table 7. Stellar parameters 



Source 




R* 




age 


E(H-K) 


Av 


Lace 






L© 


Rq 


M© 


X 10'' yr 




mag 


Lq 


Mq yr-1 


IRS2 


4.3±1.5 


2.9±0.6 


1.4±0.3'' 


0.5-r 


1.3±0.2 


20±3 


7.7±2.5 


310~^" 








1.8'' 


2" 








210-^* 


IRS5a 


1.6±0.5 


2.4±0.4 


O.SiO.l" 


0.3-0.5" 


2.9±0.2 


45±3 


~0.4 


~3 10"*" 








0.9'' 


LS*" 








~2 10"*'' 


IRS6a 


0.5±0.2 


1.5±0.3 


O.SiO.l" 


0.5-1" 


1.9±0.3 


29±5 


<0.1 


< 5 10"^ 








0.4** 


1.8'' 










HHIOO IR 


3.1±0.9 


3.6±0.7 


0.4±0.1'' 


0.1" 


1.9±0.2 


30±3 


12±2 


210"'^" 








0.75'' 


0.6' 








110"" 


IRS3 


O.SiO.l 


1.2±0.1 


0.5±0.1" 


1-5" 


0.2±0.2 


10±3 


<0.1 










0.5* 


4" 











" Masses and ages derived from the D'Antona & Mazzitelli (1997) tracks. 
'' Masses and ages derived from the Siess, Dufour & Forestini (2000) tracks. 

Table 3. HI lines in IRS2 and IRS5a Table 6. Spectral types and veiling 



Fluxes 



IRS 2 IRS5a 



^obs 


Line 


F{AF) 


F{AF) 
-icm"2 






10"^^ ergs" 


1.282 


Pa/3 


13.7(0.4) 




1.519 


Br 20 


1.6(0.5) 




1.527 


Br 19 


2.9(0.7) 




1.534 


Br 18 


3.7(1.1) 




1.544 


Br 17 


4.2(0.2) 




1.556 


Br 16 


4.1(0.4) 




1.571 


Br 15 


7.2(0.1) 




1.588 


Br 14 


6.0(0.7) 




1.612 


Br 13 


10.7(0.1) 




1.641 


Br 12 


9.3(0.6) 




1.681 


Br 11 


19.1(0.1) 




1.737 


Br 10 


41.5(0.9) 




1.818 


Br 9 


82.3(1.7) 




2.166 


Br7 


181.0(5.4) 


1.3(0.1) 



Kinomatical information 






IRS 2 


IRS5a 


A 


Line 


Vlsr 


A V 


Vlsr A V 


(/im) 




kms"''" 


kms"^ 


kms"^ kms"^ 


1.6113 


Br 13 


-13±4 


130 




1.6411 


Br 12 


-11±4 


146 




2.1661 


Br7 


-9±4 


180 


0±10 122 



The spectral types and luminosity classes, derived for 
the sources of our sample using the procedure described 
above, are given in Table 6. The effective temperatures 



Source 


ST 


Teff 




rn 


IRS2 


K2V 


4900±200 


2.9±0.5 


1.4±0.2 


IRS 5a 


K5-K7V 


4200±200 


l.OzbO.l 


O.lzbO.l 


IRS6a 


M2V 


3580±100 


O.liO.l 


O.liO.l 


HHIOO IR 


K5-M0V 


4060±300 


6.0±0.5 


3.5±0.2 


IRS3 


K5-M0 HI 


3800±100 


0.2±0.1 


O.liO.l 



Tf,ff relative to each spectral type have been derived 
adopting the conversion of spectral type into Tef / given in 
Kenyon & Hartmann (1995). We estimate that our clas- 
sification is correct inside one spectral sub-class, with the 
exception of HHIOO IR where the paucity and weakness of 
the absorption features do not allow a classification better 
than inside two spectral sub-classes. Figs I8I9I and [TOI show 
a comparison between the spectra of our sources and the 
spectra of standard stars which more closely match the 
observed spectra. Most of the sources are consistent with 
being dwarfs, as also found for other low mass YSOs (i.e. 
Greene & Lada 1997). The dwarf spectra shown by young 
stars have been recognized as an evidence for assuming 
that the observed features originate from the stellar pho- 
tosphere more than from a circumstellar disk. Disk pho- 
tospheres are indeed expected to have low gravities and 
thus should resemble giant or supergiants spectra (Greene 
& Lada 1997). 

Once defined the spectral type, the measured EWs of 
selected lines, compared with the EWs of standard spec- 
tra of the corresponding spectral type, have been used to 
derive the veiling in the K and H band. For this purpose 
we used lines not blended at the lower resolution of the 
available standard spectra, and not contaminated by pos- 
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Fig. 8. Normalized spectral segments of HHIOO IR in the 
K band compared with the spectra of a K5V (Wallace & 
Hinkle 1997) and a MOV (Greene & Lada 2001) standard 
stars. 



sible spurious atmospheric feature residuals. We also se- 
lected lines not too sensitive to the effective temperature 
to minimize the uncertainty on the spectral type deter- 
mination and also not too sensitive to the effect of stellar 
magnetic fields, like the Ti lines at 2.23/im whose EW 
can be increased by the Zeeman effect due to their high 
Lande factor (e.g. Johns-KruU & Valenti 2000). In prac- 
tice, good lines for the veiling determination resulted to 
be, in addition to the Na line at 2.209fim, when observed 
in absorption, Al at 2.1099/.tm and Ca at 2.2658/im for 
the K band, and Fe (1.6369^m), Si(1.6386Aim), in the H 
band. For IRS5a and IRS6a we do not have the medium 
resolution H band spectra, but few absorption features 
are detected in their low resolution spectra. Among them, 
Mg I at 1.7112/xm is the most suited to evaluate the H 
veiling, being not affected by strong blendings even at low 
resolution. 

The resulted veiling values are also listed in Table 6. 
From the table, it can be seen that, as expected, HHIOO 
IR and IRS2 have the largest veiling values, tk ~6 and 3 
respectively. The other considered sources show lower or 
absent veiling values. Finally, the veiling in the H band 
is about a factor of two smaller than the value. These 
veiling determinations will be discussed in Section 7.3. 



6. Stellar and circumstellar parameters 

The derivation of spectral type and veiling allows us to 
estimate other stellar and circumstellar properties from 
the available K and H photometry and source bolometric 
luminosity. The infrared color excess E{H — K) can be 
derived from the observed {H — K) color, the star color 
relative to the adopted spectral type {H — K)^ and the 
tk and ru values: 

E{H - K) = {H-K)-{H- K)o, (1) 

{H - K)o = {H- K), - 2.5 log (^:|:^) (2) 

where we have indicated with {H — K)q the dereddened 
IR color of the source. From this expression, the visual ex- 
tinction has been evaluated adopting Ay = 15.9 E{H — K) 
given by the standard interstellar reddening law (Rieke 

6 Lebofsky 1985). The derived values are listed in Table 

7 and range from Ay=l{] mag in IRS3 to Av=45 mag 
in IRS5a. Our derived values for HHIOO IR, IRS2 and 
IRS5 compare rather well with the extinction derived in 
these sources from the optical depth of the 10/im silicate 
and 3/im water absorption features, i.e. ~ 25-30 mag for 
HHIOO IR, 27 mag for IRS2 and 45 mag for IRS5(a-hb) 
(Chen & Graham 1993, Whittet et al. 1996). The low ex- 
tinction value (~10 mag) derived for IRS3 is consistent 
with values found towards field stars obscured by the R 
CrA cloud (Whittet et al. 1996). From the adopted spec- 
tral type and from the K magnitude dereddened by the 
derived Ay value, we are now able to derive the absolute 
stellar K magnitude, and thus the intrinsic stellar lumi- 
nosity: 

logU/L0 = -0.4(Mbo, - M,^i^q), (3) 
where 

Mb,i ^BC + Mk + {V- K), (4) 

In this expression, the values for the bolometric correction 
BC and for the intrinsic stellar {V — if)* colors for dwarfs 
were taken from Kenyon & Hartmann (1995). Finally, the 
stellar radius R* can be also derived given the assumed 
stellar effective temperature. Table 7 lists the L, and R* 
values for the five considered sources. Luminosities range 
between 0.3 and 4.3 Lq, and stellar radii between 1.2 and 
3.6 Rq, being therefore in the range expected for low mass 
young stars (e.g. Siess et al. 1997). Given the stellar lumi- 
nosity and its effective temperature, the sources can be po- 
sitioned in a HR diagram to have an estimate of their mass 
and age through the comparison with pre-main sequence 
evolutionary tracks. This is shown in Figure ITTl where we 
have considered the D'Antona & Mazzitelli (1997) tracks 
for stars with masses between 0.1 and 2 Mq. From these 
tracks we derive masses between 0.3 and 1 Mq and ages 
between 10^ and 10^ yr. In order to check for the sensi- 
tivity of different models used to derive the evolutionary 
tracks, we have also derived masses and ages adopting 
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Fig. 9. Normalized spectral segments of IRS2 (left) and IRS5 (right) in the K band compared with the spectrum of 
the standard star which better match the observations. The observed spectra have been rebinned at the resolution of 
the standard (~3000) for comparison. 
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Fig. 10. Normalized spectral segments of IRS6 (left) and IRS3 (right) in the K band compared with the spectrum 
of the standard star which better match the observations. The observed spectra have been rebinned at the resolution 
of the standard (^3000) for comparison. 



the Siess et al. (2000) models, which are listed in Table 
7. These last models give slightly higher masses and ages 
with respect to the D'Antona & Mazzitelli tracks, with 
mass variations which are of about 80% at worst. The 
considered standard tracks do not take into account the 
perturbation introduced by a high mass accretion rate on 
the evolutionary tracks, an effect which could be not neg- 
ligible for HHIOO IR and IRS2. Siess et al. (1997, 1999) 
have computed evolutionary tracks for stars which are still 
accreting a substantial fraction of their mass, evaluating 
the error introduced in the mass and age determination 
by taking the values expected for standard tracks. They 
show that while the determination of the (current) mass is 



not very affected by a high accretion rate, the still accret- 
ing stars are effectively younger than the standard ones, 
for ages jSl.5 10^ yr. Thus, it is hkely that HHIOG IR and 
IRS2 are actually younger than the age of ~5 lO'"^ and 10^ 
yr, respectively, that we infer from the tracks in Fig. 1111 

The knowledge of the intrinsic stellar luminosity allows 
us in principle to derive which is the fraction of excess lu- 
minosity due to accretion through a comparison with the 
source bolometric luminosity. The error on this L^cc de- 
termination is largely dominated by the precision through 
which the hhoi value is known. With regard to our consid- 
ered sources, hboi values estimated through near IR and 
IRAS photometric data are given for HHIOO IR, IRS2 
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Fig. 11. HR digram of the R CrA sources with T^ff and 
stellar luminosity derived from the analysis of the medium 
resolution IR spectra. Evolutionary tracks (short dashed 
lines) and isochrones (dot-dashed lines) from D'Antona & 
Mazzitelli (1997) are shown for stellar masses between 0.2 
and 1.7 Mq. Isochrones are reported for a - 10^ yr, b - 
5 10^ yr, c - 10^ yr, d - 5 10^ yr, e - lO"^ yr, f - 5 lO"^ yr and 
g - IQS yr 



and for IRS5 (Wilkings 1986, 1992, see Table 1). IRS5a 
is brighter and its near IR specrum is steeper than IRS5b, 
which suggests us that the bolometric luminosity of IRS5a 
dominates over the binary system: we assume a ratio 2:1 
to separate the L^oi value among the two components. For 
IRS6 and IRS3 measurements of the far IR fluxes are not 
available; however the fact that their IR spectra are not 
very steep allow us to assume that the far IR contribu- 
tion at the Ljjoi value can be negligible. We have there- 
fore given an estimate of the hhoi values for these sources 
based on their near and mid IR ground based photometry, 
and adding a correction which assumes that the spectrum 
scales as the wavelength for A larger than that correspond- 
ing to the last observed photometric flux. For IRS6 we fur- 
therly assumes that the two components a and b equally 
contribute at the total bolometric luminosity. 

From the comparison of the derived stellar luminosity 
with the Lfjoi values, it appears that only HHIOO IR and 
IRS2 have a significant excess in luminosity with respect 
to L*, which amounts to 80 % and 65 % of the Lhoi value, 
respectively. We attribute such an extra luminosity con- 
tribution to the energy released in the accretion process. 
In IRS5a, only '^20 % of the total luminosity is due to 
accretion while in the other two sources our method does 
not allow us to measure any significant accretion luminos- 
ity. Finally, the mass accretion rate can be estimated from 
the accretion luminosity. In the hypothesis that accretion 



occurs through the circumstellar disk, the total disk accre- 
tion luminosity can be written as (GuUbring et al. 1998): 



(5) 



where Ri is the inner disk truncation radius which is pre- 
dicted to be of the order of 3-5i?». The M values for HHIOO 
IR and IRS2 are ~ 210"^ and 310"^ Mq yr-\ respec- 
tively. Such values are larger than the average values de- 
rived for T Tauri stars by means of UV excess measure- 
ments (GuUbring et al. 1998), but lower than the values 
indirectly estimated in Class sources through their mass 
loss phenomena, indicating that an evolution in the mass 
accretion rate during the protostellar life is indeed very 
likely. 

7. Discussion 

The ability to separately estimate the stellar and circum- 
stellar properties of embedded YSOs, gives the possibility 
to infer in a better detail the evolutionary properties of 
these sources, to analyse their similarities with optically 
visible T Tauri stars in terms of veiling and accretion ac- 
tivity and finally to address quantitatively how their accre- 
tion properties are connected to other related phenomena 
such as jets and winds. 

7.1. Evolutionary stage of the observed sources 

The stellar parameters, accretion luminosity and mass 
accretion rates derived from the analysis of the star 
photospheric features can be used to assess the relative 
evolutionary phase between the observed sources. We 
have seen that only HHIOO IR and IRS2 have a significant 
fraction of their luminosity due to accretion, which is 
also correlated with the presence of both extreme values 
for the K and H veilings and emission line activity. On 
the other hand, IRS5a, which like HHIOO IR and IRS2 
has been classified as a Class I protostar on the basis 
of its SED between 2 and 10//m (Wilkings et al. 1997) 
does not show a significant accretion activity, as testified 
by the low Lacc/Lboi, the absence of emission lines and 
the low veiling. Presumably, this source is a T Tauri 
still embedded in its parental envelope or seen through 
a unfavourable line of sight. Finally, IRS6a and IRS3 
present characteristics of pure reddened photospheres 
where any sign of activity connected to accretion is very 
low or absent. This analysis shows that YSOs chosen to 
broadly have similar IR colors may have very different 
characteristics in terms of accretion properties. This 
finding suggests us to make an a posteriori look at the 
distribution of the IR colors of the considered sources to 
search for a finer separation between stars with different 
accretion properties. In particular the protostars HHIOO 
IR and IRS2, for which more than 50% of their luminosity 
is due to accretion, are located, in a near IR color-color 
diagram, in a region defined by colors (J-H)^4 and 
(H-K)k;2.5. In the same region can be also put the only 
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other so far identified highly accreting protostar, i.e. 

YLW15 of the p Ophiucus cloud ((J-H)>5, (H-K)==3.8, 
Greene & Lada 2002). Sources with a low or absent 
accretion activity (IRS6 and IRS3) seem to be located 
outside this box, with colors (.J-H)<4 and (H-K)<2.5. 
In this scheme IRS5a is a borderline case, which remain 
anomalous. We only note that the colors of IRS5(a+b) 
are considerably bluer than those of the other highly 
accreting protostars ((H-K)/(J-H) = 0.85 vs (H-K)/(J-H) 
~ 0.6-0.7), indicating that part of its infrared excess can 
be due to scattered radiation. In conclusion we suggest 
that the real accreting protostars can be located in a 
color-color diagram region tighter than the previously 
considered locus for Class I sources (e.g. Lada & Adams 
1992), a clue which deserves to be further investigated on 
a larger sample basis. 



It is also interesting to note how the different accre- 
tion properties of the sources correlate with their position 
on the HR diagram. HHIOO IR is the youngest source of 
the sample with an age 10''' yr, while all the others fall 
roughly on the same isochrone at t~5 10^-10^ yr. If taken 
literally, this result indicates that sources with nearly the 
same age may have very different Lacc/Lboi values, and 
thus that the age is not the only parameter regulating the 
end of the accretion process. In particular, it can be no- 
ticed the different estimated mass accretion rate of sources 
(IRS2, IRS5a and IRS6a) which appear to have a similar 
age but diff'erent masses, ranging from ~ 0.5 to ~ 1.5 Mq. 
For these sources there is an evidence that the mass accre- 
tion rate increases with the mass. Since the sources belong 
to the same star forming core, they should have started to 
collapse with the same initial cloud infall rate, which de- 
pends only on the core kinetic temperature. Therefore, it is 
not unexpected that to accumulate a diff'erent mass in the 
same ellapsing time, their mass accretion rate should have 
followed a different time evolution. Evidences for a correla- 
tion between mass accretion and stellar masses have been 
indeed reported, although with a broad scatter, for large 
samples of T Tauri and very low mass objects (Natta et 
al. 2004). 

On the other hand, we should also remind that some 

of the considered stars have a large photometric variabil- 
ity as it is indeed common in young embedded stars (e.g. 
Kaas 1999, Horrobin et al. 1997). Often variability in IR 
magnitudes is associated to variability of spectral features, 
such as HI and CO lines, whose emission is directly re- 
lated to accretion (Biscaya et al. 1997, Nisini et al. 1994). 
Therefore, if the photometric variability of the considered 
sources is related to changes of the mass accretion rates on 
short timescales, as it happens, e.g., in FU Ori-like events, 
the possibility that the three sources with the same age 
may have a similar time averaged mass accretion rate, can- 
not be excluded. 



7.2. Origin of veiling 

The observed sources present very different values of the 
IR veiling. The K veiling spans from almost in IRS3 and 
IRS6a, to a value ~6 in HHIOO IR. In T Tauri sources, rK 
values between and 2 are normally estimated (Johns- 
KruU & Valenti 2001, Folha & Emerson 1999). Higher 
values, up to 3-4, have been measured in flat-spectrum 
and Class I objects (Greene & Lada 1997, Doppmann 
et al. 2003). We however do not find a clear correlation 
with the class of the objects, at least for the three Class I 
sources of our sample which have tk values of 0.5 (IRS5a), 
3 (IRS2) and 6 (HHIOO) despite their similar spectral in- 
dex between 2 and lO/xm. For these three sources is instead 
more evident the correlation with the estimated accretion 
luminosity. In T Tauri stars, the main considered source 
of IR veiling is the emission from the accretion disk, al- 
though, given the low mass accretion rates estimated in 
T Tauri stars, disk models find some difficulties to pre- 
dict the largest observed veiling values (r^ ~ 1-2) (Folha 
& Emerson 1999, Johns-Krull & Valenti 2001). The same 
models have; also difficulties to reproduce the large IR veil- 
ing estimated in Class I sources, despite the mass accre- 
tion rate higher than in T Tauri stars (e.g. Greene & Lada 
1996). An additional source of veiling can be the emission 
from the inner circumstellar dusty envelopes. Greene & 
Lada (2002) estimate that the inner infalling envelope can 
account for about half of the K veiling they measure in the 
Class I object YLW15. Calvet et al. (1997), estimate that 
the thermal emission of envelopes with mass infall rates 
ranging from 10^"'' and 4 10^^ Mq yr^-'^, can account for a 
2.2/im veiling ranging between 3 and 9, thus in the range 
of the values we find for HHIOO IR and IRS2. The same 
models predict a H band veiling ^ 30-50% smaller than 
in the K band, consistent with our derived values. Finally, 
the small veiling observed in IRS5a, despite its significant 
IR emission excess, can be explained by assuming that its 
envelope emission is mainly due to the scattering of the 
central object emission. In this case, as pointed out by 
Calvet et al. (1997), the envelope flux contains the same 
spectral features as the object and thus does not affect the 
line EWs. Recent models for the SEDs of YSOs (Whitney 
et al. 2003) show that in Class II sources seen close to 
edge-on, extinction and scattering dominate at near-IR 
wavelengths making the SEDs of these sources similar to 
those of less evolved Class I objects. 

7.3. Correlation with accretion and ejection activity 

Our observations show that HHIOO IR and IRS2, which 
have most of their luminosity due to accretion, are also 
those exhibiting a rich emission line spectrum, with the 
strogest lines observed in the source with the highest mass 
accretion rate (HHIOO IR). 

Correlations between the accretion luminosity and the 
IR hydrogen lines luminosity have been found in T Tauri 
stars (MuzeroUe et al. 1998). Such a correlation can be 
naturally explained if the origin of HI lines is in the star 



B. Nisini et al.: Embedded sources in R CrA 



15 



magnetospheric accretion region but also if part of their 
emission come from ionized outflows, which arc ultimately 
powered by accretion. We can check if the highly accret- 
ing objects of our sample still follow the same correlation 
Lacc-Lsr-y sls T Tauri stars. Extending such a correlation 
to highly accreting embedded objects would indeed be 
important to validate the Br7 luminosity as an indirect 
quantitative tool derive the accretion luminosity of Class 
I sources. The Br7 luminosity of HHIOO IR and IRS2, 
dereddened for the extinction values listed in Table 7, are 
5.7 10~^ and 710"'* Lq, respectively. Following the fit de- 
rived by Muzerolle et al. (1998), such line luminosities 
should be consistent with accretion luminosities ~ 40 and 
3 Lhoi respectively. Considering the large errors associated 
with the Muzerolle et al. least-squares fit, our derived Lace 
values are consistent with the relationship found for the 
T Taiiri stars, extending it to about an order of magni- 
tude larger luminosities. Br7 emission is observed also in 
IRS5a with a luminosity value ~ 5 10~^ Lq. According to 
the relationship by Muzerolle et al. the corrisponding Lace 
value should be a fraction of the bolometric luminosity of 
the order of 1O~^L0, well inside the locus of the T Tauri 
star and in agreement with our roughly estimated value. 

It is usually assumed that accretion and ejection of 
matter are correlated phenomena and that this latter is 
needed to remove part of the angular momentum carried 
by the accreting material. It is therefore worthwhile to see 
if an high accretion luminosity is always associated with 
the presence of outflows. In this respect it is indicative 
that HHIOO IR is the only source of our sample having a 
clear evidence of outflowing activity, being the identified 
exciting source of the HHIOO Hebig Haro nebulosity, and 
also of the HH objects HHlOl (located in the SW) and 
HH99 (in the NE) (Hartigan & Graham 1987). A molec- 
ular bipolar outflow is also associated with this object 
as evidenced by the ammonia and HCO"*" observations of 
Anglada et al. (1989) and Anderson et al. (1997), respec- 
tively. IRS2, however, despite showing a large accretion 
luminosity value, does not show any undoubtful tracer of 
mass ejection. We have seen that the weak H2 2.12/im line 
observed in emission can be due to both the outflow and 
the disk and no other lines typical of jets, such as [Fell] 
lines, are detected in the spectrum. On a larger scale no 
optical HH object seems to be associated to this source 
(Hartigan & Graham 1987). However, the region closer to 
IRS2 has not been investigated neither in H2 imaging, to 
search for molecular jets, nor in molecular outflows radio 
tracers, so we cannot on this basis exclude that a large 
scale flow is associated with this source. 

8. Conclusion 

We have analyzed the low and medium resolution infrared 
spectra of a sample of embedded young stellar objects in 
the R GrA star forming region with the aim of constrain- 
ing their accretion properties and evolutionary stage. The 
sample includes three Class I sources (HHIOO IR, IRS2 
and IRS5) and two sources with NIR excesses in the 



Wilking et al. (1997) R GrA core survey (IRS6 and IRS3). 
The main results from this analysis are the following: 

— The low resolution spectra revealed the presence of 
strong lines in emission only in two, namely HHIOO IR 
and IRS2, out of the six observed sources. The most 
prominent detected features arc; permitted lines from 
neutral atoms, mostly from HI (Pa/3, Br7 and lines 
from the Brackett series at higher-n) but also 01, Nal 
and CI. These lines arise in a partially ionized com- 
pact and very dense gas, likely originating either in 
the star-disk accretion region, or at the base of an ex- 
panding wind. The same sources where atomic lines 
are detected also present the CO 2.3/xm bandheads 
in emission, indicating the presence of large columns 
of warm molecular gas, probably located in the inner 
sections of an accretion disk. In the medium resolu- 
tion spectra, also a weak H2 2.12/um line is detected 
in HHIOO IR and IRS2, which may be originated either 
in the circumstellar disk or in a molecular jet. 

— The medium resolution spectra show plenty of absorb- 
tion features typical of late-type stellar photospheres. 
In HHIOO IR and IRS2 such features are very weak, in- 
dicating a large amount of IR veiling. We have derived 
the effective temperature, extinction, stellar luminos- 
ity and radius of the observed sources, through a com- 
parison with standard star spectra and IR photometric 
data. This information has been in turn used to infer 
the amount of IR veiling, the accretion luminosity and 
the mass accretion rates of the objects. 

— Only in HHIOO IR and IRS2 the accretion luminos- 
ity dominates over the total bolometric luminosity 
{Lacc/Lboi ~0.8 and 0.6 respectively) and the derived 
mass accretion rates arc of the order of 3 10"^ and 
5 10~^ Mq yr~* respectively, i.e. higher by an order of 
magnitude with respect to the average values derived 
in T Tauri stars. In contrast, in IRS5a only ^2% of 
the luminosity is due to accretion. In general, we found 
that there is a correlation between the accretion lumi- 
nosity, the IR veiling and the emission line activity of 
the sources. 

— If compared with standard evolutionary tracks for pre- 

main sequence evolution, the derived stellar luminosi- 
ties and effective temperatures indicate that HHIOO 
IR is the youngest of the sources, with an age ~10^ yr. 
On the other hand, IRS2, IRS5a and IRS6a have about 
the same age (~5 10^- 10^ yr) despite the large differ- 
ences in their accretion properties. For these sources 
there is a hint that the mass accretion rate, at a given 
age, may depends on the accumulated stellar mass, a 
result not totally unexpected if one assumes that the 
mass accretion rate is a function of the mass of the 
envelope. We cannot however exclude that variations 
on short timescales of the accretion rates, related to 
the sources infrared variability, may be responsible of 
the derived differences in the M values. 

Our analysis has shown that sources broadly defined as 
Glass I protostars, may indeed be very different in terms of 
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accretion activity. This indicates how the definition of the 
evolutionary stage of deeply embedded YSOs by means of 
IR colors needs to be more carefully refined. In particular, 
on the basis of our results we suggest that the highly ac- 
creting protostars (i.e. those for which Lacc/Lboi > 50%) 
can be located in a color-color diagram region character- 
ized by near IR colors (J-H)^4 and (H-K)^2.5, a clue 
which remains to be further investigated on a larger sam- 
ple basis. Low resolution IR spectroscopy can give a better 
diagnostic about the accretion activity of young embed- 
ded sources, evidenced by the presence of emission lines, 
HI recombination lines and CO overtone emission in par- 
ticular. However, we have shown that the combination of 
low and high resolution spectroscopy may be effective in 
constraining both the stellar and circumstellar properties 
of even highly veiled objects, allowing to directly derive 
the accretion luminosity and mass accretion rates back to 
the first 10~^ yr of protostellar evolution. 

Appendix A: IR features selected for spectral 
classification 

In the observed region of the K band we have identified 
two spectral segments, around 2.12 and 2.23/im, where 
several features can be adopted for spectral classification. 
This can be seen in Figs lA.ll and IA.2I where standard 
spectra of K and M spectral types, taken from the Wallace 
& Hinkle (1997) catalogue, are plotted. In the spectra of 
dwarfs, around 2.12/im Si at 2.1360 /imand Mg at 2.1066 
fim weaken going from early K to M types and hence their 
ratios with respect to the Al lines (2.1099, 2.1169/im), 
whose EWs remain fairly constant, are sensitive to the 
effective temperature. 

However, in the giant spectra such ratios are less sensi- 
tive to the effective temperature, and therefore a K5V star 
can be easily confused with a M2III photosphere. In order 
to distinguish also the spectral class, we can use other lines 
around 2.23^m. For example the ratios Ti(2.2217, 2.2240 
/xm)/Fe(2.2266 ^lm) and Sc(2.2401/im)/Fe(2.2386Aim) sig- 
nificantly change between dwarfs and giants in K stars 
(Fig. US). 

For the determination of the spectral type and lumi- 
nosity class in our source sample we mainly adopted the 
above EW ratios in the K band, and checked a posteri- 
ori that the other observed features were well reproduced. 
The use of the features observed in the medium resolution 
H band spectral segment for the spectral type determina- 
tion is limited due to the strong blending affecting most of 
the lines, and owing to the fact that the strongest and less 
blended lines, like Fe I at 1.637/im and Si I at 1.639/zm, 
are not very sensitive to the effective temperature. On the 
other hand, this fact makes them very useful to determine 
the value of the veiling in the H band once the source 
spectral type have been derived from the K band features 
analysis. 
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Fig. A.l. Spectra of dwarfs (on the left) and giant (on 
the right) standard stars taken from Wallace & Hinkle 
(1997) in a spectral segment around 2.1/im, which includes 
spectral features of Mg, Al and Si. 



Livia Origlia for useful discussions on spectral classification and 
Tom Geballe for helpful comments on the unidentified features. 
This research has made use of NASA's Astrophysics Data 



I I I I I I I 



I I I I I I I 



I I I I I I I I I I I I I I 

Ti Ti FeTi Ti 
III 



Fe Sc" 




I ' ' ' ' I ' ' ' ' I ' ' ' ' I 



I I I I I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I I I 

rr\ Ti FeTi Ti Fe Sc " 

I ..^1 . nl J, 



2.22 2.225 2.23 2.235 2.24 



2.22 2.225 2.23 2.235 2.24 
A(/im) 
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(1997) in a spectral segment around 2.3/im, which include 
spectral features of Ti, Fe and Sc. 
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